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We present the results of three-dimensional molecular dynamics simulations of vortices which indicate that, 
for B larger than the matching field, the enhanced pinning effectiveness of splayed columnar defects relative to 
vertical columnar defects can be explained in terms of the existence or absence of channels through which the 
vortices can flow without encountering defects. 



1. Introduction 

In recent years, experiments have shown that 
the critical currents of high-temperature super- 
conductors can be greatly enhanced by the intro- 
duction of long, straight defects into the mate- 
rial, particularly if these defects are splayed at 
some small angle with respect to the crystalline 
c-axis|Q]. These experiments have generally been 
performed such that the magnetic field is less 
than the matching field, although recent efforts 
have explored the opposite case^,||. Explana- 
tions of this effect are incomplete. It has been 
suggested that random misalignment of the pins 
inhibits large-scale, low-energy excitations, and 
that entanglement further prohibits vortex mo- 
tion P]. At the same time tilted columns yield an 
increase in the bending energy of the vortices and 
promote hopping at defect crossings. Our results 
indicate that the physics of splay enhancement 
may be missing an important ingredient. We find 
that splaying columnar defects slightly closes off 
channels through which vortices can flow, produc- 
ing enhanced critical currents. 

2. Model 

We have conducted molecular dynamics sim- 
ulations of driven, interacting vortices in three 
spatial dimensions. The forces on the vortices 
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Figure 1. Resistivity p as a function of driv- 
ing current density J for point (circles), col- 
umn(squares), and 3.6° splay (diamonds). 

result from the vortex-vortex interaction, the 
vortex-defect interaction (modeled as a short- 
range attractive potential), elastic bending, ther- 
mal Langevin forces, damping, and an exter- 
nal Lorentz force. The energy scales are set by 
V VO rtex = ( < I > o/47rA) 2 , with $o the flux quantum 
and A the magnetic penetration depth. Periodic 
boundary conditions are imposed in the a-b plane 
so as to maintain a constant flux density, while 
open boundary conditions are imposed along the 
c-axis. We have chosen to simulate 30 vortices on 
20 planes in a 16A x 16A periodic cell, with 20 
defects per plane. At this vortex density, the vor- 
tex repulsion is strong enough to prohibit vortex 
entanglement. The defects were arranged cither 
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Figure 2. J C (T) for point (circles), column 
(squares), and 3.6° splay (diamonds). J c — at T 
= 0.225, 0.250, and 0.325 for point, column, and 
splay, respectively. The error bars are roughly 
twice the point size. 

as vertical columns, columns splayed at an angle 
O with respect to the c-axis, and uncorrelated 
point defects. The defect radius is taken to be A, 
the inter-planar spacing is 12A, and the coherence 
length £ is 24 A. 

3. Results 

A typical result is shown in Fig. 1 for columnar, 
point, and splayed defects at low temperature. 
The current density is normalized by the BCS 
depairing current density, while the resistivity is 
normalized by the Bardeen-Stephen flux-flow re- 
sistivity. For columnar defects at very low cur- 
rents, vortices are either pinned to defects or held 
in place by their mutual repulsion. As the current 
density increases, some of the interstitial vortices 
shear off into channels of vortex flow. Point de- 
fects leave no open channels for vortices, but the 
scattering of the defects reduces their ability to 
pin vortices. Splayed defects also leave no chan- 
nels, and, at small angles, are able to effectively 
pin entire vortices. 

Fig. 2 shows the critical current density J C {T) 
for each configuration, where J c is defined by a 
threshold criteria p{J c ) = 0.05. For point de- 
fects, J c drops rapidly with increasing tempera- 
ture, reflecting their weak pinning efficiency. For 
columnar defects, J c drops slowly with increasing 
temperature as more channels open for the vor- 



Figure 3. Critical current density as a function of 
splaying angle from the c-axis for T = 0.025. 

tices. Splayed defects have the largest J c for all 
temperatures. 

We consider the variation in J c for different 
splay angles at low temperature in Fig. 3. J c 
shows an enhancement at small angles due to the 
reduction of available channels. For large angles, 
the vortices can no longer accommodate to the de- 
fects, and splayed defects become similar to point 
defects. The value of the maximum, roughly 5°, 
compares extremely well to the number observed 
in measurements in YBa2Cu307j^]. 

In summary, our simulations show that the crit- 
ical current enhancement from splay is the result 
of reduced vortex channel flow rather than entan- 
glement for this vortex density. We are exploring 
other regions of parameter space to determine op- 
timal conditions for splay enhancement. 
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